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In numerical studies of turbulence, hyperviscosity is often used as a tool to extend the inertial 
subrange and to reduce the dissipative subrange. By analogy, hyperdiffusivity (or hyperresistivity) 
is sometimes used in magnetohydrodynamics. The underlying assumption is that only the small 
scales are affected by this manipulation. In the present paper, possible side effects on the evolution 
of the large scale magnetic field are investigated. It is found that for turbulent flows with helicity, 
hyperdiffusivity causes the dynamo-generated magnetic field to saturate at a higher level than normal 
diffusivity. This result is successfully interpreted in terms of magnetic helicity conservation, which 
also predicts that full saturation is only reached after a time comparable to the large scale magnetic 
(hyper)difi'usion time. 

PACS numbers: 47.11. -(-j, 47.27.Ak, 47.65.-fa, 52.65.Kj 



In theoretical studies of Navier-Stokes turbulence the 
ordinary viscosity operator, i^V^it, is sometimes replaced 
by (— 1)"~^z/„V^"m, where i^n is a hyperviscosity of or- 
der n. The use of hyperviscosity has the advantage of 
making the transition from the inertial subrange to the 
viscous subrange shorter At the same time, however, 
it has the notorious disadvantage of tempering possibly 
major parts of the inertial subrange. An example is the 
so-called bottleneck effect that leads to significantly shal- 
lower power spectra at high wavenumbers near and before 
the viscous subrange Q. When early simulations using 
the piecewise parabolic method showed such bottleneck 
effect it was unclear whether this effect was real or 
just a consequence of hyperviscosity. In the context of 
magnetohydrodynamics, recent comparisons between di- 
rect and hyperviscous simulations point now to the latter 
possibility 

The use of hyperviscosity is indeed quite popular in 
studies of hydromagnetic turbulence where, in addition 
to viscosity, the ordinary magnetic diffusivity is replaced 
by hyperdiffusivity [§, ||, ^ 0. A strong artificial bottle- 
neck effect occurs when hyperdiffusivity is used Q . This 
is particularly clear in two dimensions at very high reso- 
lution [^, although a weak bottleneck effect occurs even 
without any hyperdiffusive effects |^ . The use of hyper- 
viscosity and hyperdiffusivity has also led to significant 
controversy Q in models of the Earth's dynamo ^ . At 
the center of the controversy is the effect of hyperviscos- 
ity on the asymptotic behavior at small Ekman number 
(low viscosity). This can also affect conclusions regard- 
ing the relative importance of the Lorentz force, and the 
relevance of Taylor's constraint, both matters of great 
importance for geodynamo theory. There are also exam- 
ples where the use of hyperdiffusivity has moved dynamos 
from an a^-regime toward an ail-regime [ pO| . 

Generally speaking, hyperviscosity and hyperdiffusiv- 
ity can lead to rather ill-understood behavior that tends 
to diminish its potential advantages. It is therefore im- 



portant to clarify exactly how hydromagnetic dynamos 
are affected by this approach. Here, we concentrate on 
the effects of (magnetic) hyperdiffusivity. This is done in 
the context of MHD turbulence; future work will look at 
this effect on geodynamo models. 

In the present paper we show that, if the fluid motions 
are helical, hyperdiffusivity can lead to artificially en- 
hanced saturation amplitudes of the nonlinear dynamo. 
It is at first glance somewhat counterintuitive that the 
properties of the large scale field should depend on the 
details of the diffusion operator, which is supposed to af- 
fect only the small scales. In recent years, however, there 
has been mounting evidence that large scale dynamos, 
which usually involve helicity, do depend on the micro- 
scopic diffusion [|ll| . This property is related to magnetic 
helicity conservation, which permits magnetic helicity to 
change only on a resistive time scale; see Ref. here- 
after referred to as B2001. These processes can be seri- 
ously affected by the use of hyperdiffusivity. Obtaining a 
detailed understanding of the associated artifacts is cru- 
cial before hyperdiffusivity can be taken as a useful tool 
in dynamo simulations. We emphasize that the sensitiv- 
ity to the use of hyperdiffusion reported in the present 
paper is peculiar to large scale dynamos and does not 
apply to small scale dynamos. 

The magnetic field evolution is governed by the induc- 
tion equation. 



dB 



= V X (m X B) + (-1)"-177„V2"B, 



(1) 



with V • i? = 0, rjn = const, and n = 1 for ordinary 
magnetic diffusivity. This equation has to be integrated 
together with the momentum and continuity equations 
which are, for an isothermal compressible gas with con- 
stant speed of sound, Cg, 



Dm 



— = -c^Vlnp 



J X B 



(2) 
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FIG. 1: Evolution of magnetic energy spectra in Run A (with 
ordinary magnetic diffusion) in equidistant time intervals at 
t = (lowest curve), t = 80 (next one higher up), until t = 
1200 (peaking at the very top left). The dotted lines give the 
kinetic energy. 



where D/Dt — d/dt + u ■ is the advective deriva- 
tive, Fvisc = p(V^M + |W • u) is the viscous force, 
fj, = const is the dynamical viscosity, Cg is the isothermal 
sound speed, J = V x is the current density, and 

/zq is the vacuum permeability. As in B2001, the forc- 
ing function f(x,t) is a randomly chosen polarized wave 
taken from a band of wavenumbers around the forcing 
wavenumber k{. The direction and phase of / change 
randomly at each time step. We solve the equations in 
a periodic domain of size L^, where L — 2tt. We adopt 
nondimensional units by measuring density in units of the 
average value, po = {p)^ length in units of l/fci, where 
/ci = 27r/L = 1 is the smallest wavenumber, velocity in 
units of Cs, magnetic field in units of Cg-^poMOi and time 
in units of (csfci)"^. 

For numerical solutions, the smallest possible value 
of rjn is given by the condition that the mesh magnetic 
Reynolds number, based on the smallest scales resolved, 

R^r^''^ = ^rms/(?7„4r') 

is of order unity. Here, /cNy — t^N/L is the Nyquist 
wavenumber of a mesh with TV points. 

As a reference model, we consider a calculation with 
ordinary magnetic diffusivity, ryi = 10~^ and a forcing 
wavenumber fcf = 27 (Run A). We adopt a dynamical 
viscosity of /i = 10^^, so the magnetic Prandtl num- 
ber is 100. In Fig. |l| we show spectra at different times. 
Consistent with earlier results (B2001), magnetic energy 
grows owing to dynamo action with spectral peaks at 
the forcing scale, fc = fcf , and at some intermediate scale. 



fc « 9. By t « 480 (seventh curve from the bottom of 
Fig. |l|) the intermediate scale field has reached equiparti- 
tion with the kinetic energy; the field continues to grow, 
however, and now evolves towards larger scales under a 
envelope. 

Qualitatively similar behavior is found for the hyper- 
diffusive Run B, with 772 = 3 x 10^®; see Fig. || This case, 
which has similar diffusion at the Nyquist wavenumber 
to Run A, also exhibits a secondary peak at some in- 
termediate wavenumber. Following B2001, cf. Eq. (36)- 
(39) therein, we identify this with the wavenumber where 
growth rate of a corresponding dynamo model is max- 
imum. In the initial kinematic stage, the position of 
the secondary peak is approximately constant in time 
(fcmax ~ 9 for Run A and fcmax ~ 14 for Run B). When 
the magnetic energy reaches equipartition with the ki- 
netic energy, Skin, the secondary peak begins to travel 
toward smaller fc. A reasonable fit to this migration is 
given by 

^max = Q!tiav(i " ^sat), (5) 

where the parameter atrav characterizes the speed at 
which the intermediate peak travels. The values obtained 
from fits to these runs are listed in Table |, together 
with some other parameters defined below. Note that 
for Run C, with a lower hyperdiffusivity 772 ~ 10^^, the 
speed of the peak has decreased even further. This sug- 
gests that even at intermediate scales (i.e. fc less than the 
initial value of fcmax), the dynamo process is resistively 
limited. 

TABLE I: Summary of the runs discussed in this paper. Note 
that Qtrav decreases with decreasing values of A is the 
kinematic growth rate of the magnetic energy, H and M are, 
respectively, magnetic helicity and energy during the kine- 
matic stage, and the parameter ^skin (defined below) gives an 
approximate upper bound for In = \H\/{2^oM). 

Run N n 1]„ A fcf Qtrav £h ^skin 

'A 120-* 1 0.047 27 1.1 x 10"'' 0.035 0.065 

B 120^ 2 3 X 10"* 0.070 27 7.3 x 10"'' 0.018 0.025 

C 120^ 2 10"** 0.082 27 3.6 x 10"* 0.005 0.013 

D 30^ 2 10"* 0.078 3 - 0.08 0.15 



Since it is difficult to evolve a simulation at this reso- 
lution over a full large scale diffusion time, ~ (77„fc^")~^, 
we now compare with a low resolution run with only 30'^ 
mesh points, 77„ = 10""* (Run D). The parameters of this 
run are chosen so that the factor fc^""^, which appears 
in the theory below, is consistent with Run A. The large 
scale magnetic field shows a very prolonged saturation 
phase after the saturation of the small scale field (and 
the equipartition of magnetic and kinetic energy), finally 
equilibrating only after approximately one large scale dif- 
fusion time; see Fig. ^. 

As in B2001, we can interpret the slow saturation be- 
havior in terms of the magnetic helicity equation. We 
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FIG. 2: Evolution of magnetic energy spectra in the hyper- 
diffusive Run B in equidistant time intervals at f = (lowest 
curve), t = AQ (next one higher up), until t — 2080 (peaking 
at the very top left). The dotted lines give the kinetic energy. 




FIG. 3: Evolution of large scale magnetic energy, normalized 
by the kinetic energy, for Runs A, B and D. The dotted lines 
give the result expected from Eq. (p^, as is explained below; 
these lines show strong fluctuations because we have used here 
the actual value of M{. For Run D, Mf is also shown (grey 
line at the very bottom). The inset shows the early kinematic 
evolution for Runs A-D in a semi-logarithmic representation. 



begin with the uncurled induction equation, 
dA 



dt 



uxB + (-l)""^r7„V^"A - V 



(6) 



where A is the magnetic vector potential, with B = 
V X A, and is the electrostatic potential which can 
be chosen arbitrarily without affecting B. For a periodic 
domain of volume V, the magnetic helicity. 



H 



A- BdV = {A- B)V, 



(7) 



where angular brackets denote volume averages, is inde- 
pendent of the choice of 0. (For the simulations we take 
(/> = 0.) Dotting Eq. (|) with A, and Eq. (|) with B, 
adding the two and averaging, yields 



-l)"2ry„((v2"A).iJ) 



(8) 



Surface terms are absent because of periodic boundaries. 
For this reason, and because V • B = 0, the right hand 
side of Eq. (H) becomes —2rji{J ■ B) when n — \. 

We now proceed analogously to B2001. Firstly, in 
the steady state, {A ■ B) must be constant and so 
((V^"A) • B) must vanish. This happens in such a way 
that there are contributions from the forcing scale and 
the large scale such that the two terms have opposite sign 
and cancel. We calculate the large scale field, B = V x A, 
by Fourier filtering around k = \ (using integer bins) , and 
the small scale field, b = 'Vx a, as b = B — B. The char- 
acteristic wavenumbers of these scales are k{ (for forcing 
or fluctuating scale) and ki (for smallest wave number). 
Making use of the fact that the magnetic field is nearly 
fully helical at small and large scales, we have 



(-1)" ((V2"a) • b) « ±k^"'\b^) 



(-1)"((V2"A) -B) w 



(9) 



(10) 



where the upper and lower signs are for positive and neg- 
ative signs of the kinetic helicity of the forcing, respec- 
tively. Thus, the ratio of large to small scale magnetic 
energies, i.e. the degree of superequipartition, is 



Ml 



(62) 



^ >1. 



(11) 



For fcf — 3, normal diffusion {n = 1) gives superequipar- 
tition by a factor of 3; see B2001. For Run D {n — 2) 
we should have superequipartition by a factor of 27. The 
numerical result (Fig. ||, where M{/{E\^in) ~ 0.5) gives 
Ml /Mi « 44. As we explain below, this has to do with 
the fact that the estimates in (||) and (|l^) are not quite 
accurate. Nevertheless, the effect of hyperdiffusivity on 
the level of superequipartition is clear. 

Analogously to B2001 we can also calculate the asymp- 
totic saturation behavior by using Eqs @j and (^ and 
assuming equipartition at small scales, (b^) « {pu'^), 
which is expected to hold after the time tg when the small 
scale field has saturated. (The fact that these solutions 
satisfy M{/ (-Ekin) ~ 0.5, rather than strict equipartition, 
does not affect the following.) This gives 



which has the solution 



2n-l 



(#) = (b') 



2ji-l 



B 



1 - e 



2rink 



2"-l(62\ 



-2r;„fef*(t-ta) 



(12) 



(13) 
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In the early saturation phase, we have 



(14) 



Thus, Run A (r/i = 10"-* and fcf = 27) and Run D (772 = 
10^"* and kf — 3) should exhibit the same saturation 
behavior; this can be approximately verified from the 
early saturation behavior visible in Fig. ^. 

The magnetic field in Run D actually saturates some- 
what faster than suggested by Eq. (p^. Again, this is 
explained by the observation that the estimates for the ef- 
fective values of the wavenumbers in (|^) and (|lO|) are not 
accurate. Good agreement can be achieved if, instead, 
we use Eqs @ and to calculate effective wavenum- 
bers, fci —f fci^off and fcf — > kf^cSi for the large and small 
scale fields, respectively, and if we use the actual values 
for the small scale magnetic energy. Aft (which show a 
long-term trend, but is also fluctuating). For Run D we 
find fci^cff = 1-3 and /cf^off — 4.6. Such an enhancement 
results from hyperdiffusivity which increases the relative 
contributions from higher harmonics. The modified ver- 
sion of Eq. (nM is then 



Ml = Mf 



tf.off 



ki_, 



2n-l 



_ p-2i)„/cifei"f/ (t-ts) 



(15) 



The evolution predicted by Eq. (|T^) is shown as dot- 
ted lines in Fig. ||. Note that the time taken for satu- 
ration is dependent upon the large scale hyperdiffusion 
time (77„fcJ"jj)~^, but that since the large scale field is of 
approximately unit wavenumber, fci_oft ~ 1, the hyper- 
diffusivity has very little effect, decreasing this time only 



slightly (cf. the true large scale diffusion time for this 
value of 77). In this respect, hyperdiffusivity is behaving 
exactly as we would wish; allowing us to attain low 77 
at lesser computational expense, and with little effect on 
the physical behavior. 

We note that during the kinematic phase the dynamo 
is still growing on a fast dynamical time scale. At this 
stage, the net magnetic helicity remains close to zero, 
as it must for the high magnetic Reynolds numbers un- 
der consideration. Berger's inequality fl^ gives an up- 
per limit for the growth of magnetic helicity, derived by 
bounding the right hand side of Eq. (^) via the square 
root of Joule dissipation and magnetic energy. In the 
presence of hyperdiffusivity this inequality is 



\H\/{2^loM) < at 



skin; 



(16) 



where 4kin = (277nfcf "^^/A)^/^ is a modified skin depth, 
A is the kinematic growth rate of the magnetic energy, 
and a is a coefficient of order unity. From Table | we 
see that this constraint is indeed well satisfied during the 
kinematic growth phase. 

The present results have demonstrated that hyperdif- 
fusivity can have profound effects on dynamos with he- 
licity. The modifying effects are well understood, which 
makes the use of hyperdiffusivity an efficient tool for nu- 
merical studies. This has allowed us here to show that 
helical dynamos saturate resistively both on large and 
intermediate scales, but not on small scales. 

Use of the PPARC supported supercomputers in St 
Andrews and Leicester (UKAFF) is acknowledged. 
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